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We analyze a possible extension of Gioia and Chakraborty's momentum transfer model of friction 
in steady turbulent pipe flows (Phys. Rev. Lett. 96, 044502 (2006)) to the case of time and/or space 
dependent turbulent flows. The end result is an expression for the stress at the wall as the sum of an 
steady and a dynamic component. The steady part is obtained by using the instantaneous velocity 
in the expression for the stress at the wall of a stationary flow. The unsteady part is a weighted 
average over the history of the flow acceleration, with a weighting function similar to that proposed 
by Vardy and Brown (Journal of Sound and Vibration 259, 1011 (2003); ibid. 270, 233 (2004)), but 
naturally including the effect of spatial derivatives of the mean flow, as in the Brunone model (B. 
Brunone et al., J. of Water Resources Planning and Management 236 (2000)). 



1. INTRODUCTION 

When a fluid flows through a pipe of circular section and radius R it experiences a pressure drop per 
unit length of pipe dp/dx = —2t/R, vifhere t is the stress at the wall, r has units of energy density. For 
time-independent flows it is commonly parameterized as in the Darcy-Weisbach formula 

r^^pU' (1) 

where p is the density of the fluid and U the average flow velocity across a given transverse section of the 
pipe. The coefficient / in cq. ([ij is the so-called friction factor [l|-|5|. 
For a given pipe, the friction factor is a function of Reynolds number 

Re.^ (2) 

V 

where v is the kinematic viscosity of the fluid (as distinct from the dynamic viscosity p = pv). f presents 
three power law-like regimes separated by transition regions. For laminar flows (Re < 10'^) we have the 
Poiseuille Law 

/==64Rc"^ (3) 
for developed turbulent flows (10"^ < Re < 10®) it obeys the Blasius Law 



/ = 0.3164 Re-i/" (4) 

and for larger values of Reynolds number it converges to an asymptotic value determined by the pipe 
roughness e according to the Strickler Law 



1calzetta@df.uba.arl 



2 



/ = 0.142 6^/3 (5) 

e is defined as the characteristic ratio of waU protrusion height to pipe radius. 

G. Gioia and P. Chakraborty @ have presented a theoretical model whereby the three power law regimes 
of the steady flow friction factor are easily derived from the Kolmogorov spectrum of homogeneous, isotropic 
turbulence. See also We shall refer to this as the momentum transfer model (MTM). The original 

presentation of the MTM made contact with the Blasius and Strickler (for rough pipes) asymptotic regimes, 
while the transition from the Blasius to the Poiseuille regimes @, was discussed in [l5[ . The connection 
between wall friction and turbulent spectrum can be used in two directions. If the turbulent spectrum can 
be found, then it can be used to derive properties of the wall friction. This has been the spirit of the original 
presentation In [T]| this approach has been used to explore friction in two dimensional flows, and in [TBj 
to investigate drag reduction by polymer additives. 

Alternatively, we may use properties of friction to derive certain features of the turbulent spectrum. For 
example, in jlOf the scaling properties of the friction coefficient are used to investigate intcrmittency in the 
turbulent part of the flow. 

In this paper, we aim to combine both approaches to study friction in time-dependent and non- 
homogeneous turbulent pipe flow. We shall show that a relatively simple generalization of the MTM yields 
a reasonable expression for the unsteady friction in the somewhat academic case of inhomogeneous, time- 
independent flows [l^. However, the same generalization fails to describe correctly unsteady friction for 
time-dependent, homogeneous flows. We shall consider on a phenomenological basis how the MTM should 
be modifled to capture at least some basic elements of unsteady friction in this second case. If the philosophy 
of the MTM is upheld, this sheds light on the way a time-dependent mean flow excites turbulent fluctuations. 

The subject matter of time dependent turbulent pipe flows has considerable interest, both theoretical and 
practical |17H2l| . The basic framework of analysis depends upon the separation of the flow velocity and 
pressure into a mean component and a turbulent fluctuation. The action of the turbulent fluctuations on the 
smooth part of the flow is encoded into the stress at the wall by a suitable constitutive relation. To make 
things simpler, in most cases of interest the turbulent part of the flow may be regarded as incompressible, 
since characteristic velocities are much lower than the speed of sound. This allows one to bring the extensive 
lore on friction in incompressible pipe flow to bear on the problem. We shall also assume that the mean part 
of the flow is parallel. For an incompressible flow this means the mean flow is independent of position along 
the pipe. Wc shall relax the incomprcssibility assumption when discussing inhomogeneous flows. 

Another important simpliflcation is available at very high Reynolds number Re. In this case, mean flow 
may be regarded as essentially one dimensional, and is determined by the equation of state for the fluid and 
the conservation laws for mass and linear momentum, namely 



djpU) d{pU^) ^_dp_2T 
at dz dz R ^ ' 

where r is the stress at the wall. It is convenient to build in the known friction factor for steady flow by 
writing t = + Tf. Here r, is the steady part of the friction factor given by eqs. is the dynamic 

contribution to friction. 

Most work on time - dependent turbulent pipe flows simply neglect Tt p2j . This approach is sensible at 
low frequencies, but it is bound to fail eventually [23l425| . Indeed that is the case in time-dependent laminar 
flows, where the time dependent friction may be computed with great accuracy [1^. 

In the literature |27H32| |. it is common to group the different approaches to friction in time-dependent, 
inhomogeneous flows under two broad headings, namely the convolution method (CM) [2^, [33l - l35| and the 
instantaneous acceleration method (lAM) [36|, which in turn may be derived form extended irreversible 
thermodynamics (EIT) [13 • These approaches are described in next Section. 
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To apply the MTM under unsteady conditions we shall follow the same strategy we have used successfully 
to study drag reduction from polymer additives in |15| . Namely, we shall substitute the turbulent dynamics 
in the stationary case by a linear stochastic equation which reproduces the Kolmogorov spectrum. Then we 
shall write down a generalization of this equation to include a nontrivial background velocity. We shall solve 
the generalized Langevin equation to compute the spectrum under the new conditions, and thereby derive 
the friction by following the MTM prescription. 

The fact that the correlations of turbulent fluctuations may be cast as stochastic averages over a suitable 
Langevin dynamics is just a particular application of a basic theorem in nonequilibrium field theory (38l - l40j . 
The problem is that to really derive this stochastic dynamics from the Schwinger-Dyson equations we need 
a detailed knowledge of the fluctuation self correlations that is presently unavailable, even for steady flows. 
Therefore the best we can do is to motivate a particular ansatz for the Langevin equation, by using the few 
known facts about the spectrum, physical common sense and Occam's razor. 

As a matter of fact, this approach works well for time- independent, inhomogeneous flows. Generally 
speaking, the MTM understands wall friction in steady flows as the incoherent sum of contributions from 
turbulent eddies of different sizes. The contribution from each eddy is proportional to its characteristic 
velocity, which in turn is proportional to the mean flow velocity. A natural generalization of this scenario 
to inhomogeneous flows is to assume that the eddies contributing to friction at a given position have been 
created at different points upstream, and then advectcd by the mean flow (being slowly dissipated in the 
process) [4l| - |43j . Following this line of argument, it is natural to assume that the relevant eddy velocity is 
proportional to the mean velocity at the point of creation, rather than at the point where wall friction is 
being computed. With just this set up we shall be able to arrive at a Brunonc-likc formula for the unsteady 
friction. 

However, the same approach fails to reproduce the Vardy and Brown CM model of turbulent friction in 
fully rough flows, which we shall take as paradigmatic, in the case where the mean flow is homogeneous 
but time dependent. Consider for example the case of an uniformly accelerating flow [il], and let us call 
"now" the time when we want to compute wall friction. At any time in the past mean velocities are lower 
than now, and therefore the wall friction computed as above will be necessarily lower than the friction for a 
steady flow with the present instantaneous mean velocity. Vardy and Brown's CM, and most of the models 
in the literature, predict just the opposite. Confronted with this fact we shall seek a minimal modiflcation 
of the Langevin equation which reproduces at least the basic elements of the CM, without including more 
free parameters than already present in the Vardy and Brown approach. Although this will not teach us 
much about unsteady friction, since it is essentially a transcription of the Vardy and Brown model into the 
Gioia and Chakraborty language, it is potentially revealing concerning the way unsteady pipe flows excite 
turbulence. 

The rest of the paper is organized as follows. In next section we include some necessary background 
material. We first review the basic approaches to unsteady friction [2^, [23| , with emphasis on those of Zielke 
[26| . Vardy and Brown [sl, and Brunone et al. [1^. Then we introduce the MTM, following closely the 
original presentation of [g. We finally describe the stochastic approach to the MTM previously introduced 
in [1^. The following two sections are the bulk of the paper, since there we generalize the MTM to time- 
independent inhomogeneous and homogeneous, time-dependent situations, respectively. We conclude with 
a few final remarks. 



2. BASIC BACKGROUND 



2.1. Main approaches to transient friction 



In this subsection we shall briefly review the main approaches to transient friction in the literature. 

Zielke [1^ provides a thorough discussion of transient friction in laminar flows. By assuming the flow 
may be regarded as parallel and incompressible, the full set of continuity and Navier-Stokes equations [U 
is reduced to a single equation for a radially dependent mean velocity u [z, r, t], related to U [z, t] in eqs. © 
and ([7]) through 
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U[z,t] = — / rdru[z,r,t] (8) 



If we further neglect z-derivatives and nonlinear terms, this single equation reads 



du v 1 d du 1 dp 
dt r r dr dr p dz 

Average of this equation over a cross section yields cq. ([7]) with the identification 

du 

T^-vp—[z,R,t] (10) 
Because equation ^ is linear and z appears only as a parameter, it admits a solution where 

u [z, r,t]^- f dt' g[r,t~ t'] ^ [z, t'] (11) 

P J-oo oz 

Averaging this equation over a cross section yields a similar relationship 

U[z,t] = - f dt' G[t~t']^[z,t'] (12) 

P J -oo 

Wc may consider this as an integral equation for p for a given U , invert it and rewrite eq. ()11|) as 

u[z,r,t]= [ dt' gi[r,t-t']U[z,t'] (13) 



and finally use this equation and (jlOp to obtain a linear relationship between the stress at the wall and the 
mean velocity. Thus unsteady friction is naturally expressed as a convolution over the history of the mean 
flow. 



Tt[z,t]^pJ dt'W[t^t']-^U[z,t'] 



(14) 



The CM attempts to retain this pattern in the turbulent case, which requires some ad hoc assumptions pSj . 
For example, the Vardy and Brown weighting function is derived by writing a linear equation for the mean 
velocity along the lines of eq. (|9]), but with an effective viscosity that depends on the radial distance to the 
pipe. Although more complex, the analysis follows the same steps as Zielke's. The main result of the analysis 
is the "weighting function" W [t] . The particular weighting function derived by Zielke has been found to 
be an accurate description of transient friction both for laminar and low Reynolds number turbulent flows 
[23} . Vardy and Brown [s^ [s^l and others [i^, \^ have derived weighting functions appropriate for fully 
developed turbulent flows. Most proposals share many features, like a pronounced peak or else an integrable 
singularity as t — > and an exponential fall off for large t. We shall take the Vardy and Brown weighting 
function as paradigmatic [33 |. The Vardy-Brown weighting function is best introduced through its Laplace 
transform 



(15) 
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which in time domain becomes 



WvB M^^J^e-^"* 



(16) 



where w is the viscous dissipation frequency 



UJ 



i?2 



(17) 



The constants A and C depend on Reynolds number and roughness. For fully rough flow they are given by 



'g \ 0.39 

A = 0.0206 VRe ( - 



(18) 



C = 0.352 Re ( - 



(19) 



Observe that Re in eqs. (|18p and ((T^ is not the instantaneous Reynolds number eq. 0, but rather a 
fiducial Reynolds number, e. g., the Reynolds number of a preexisting steady flow [s^l. By construction, Re 
is space and time independent. 

While the CM represents unsteady friction as an average over the history of the mean flow at a flxed point, 
the lAM builds unsteady friction from time and space derivatives of the mean flow at a given time. The 
Brunone ansatz for the transient friction reads 13611 



Tt [z,t] = kspR 



dU 
'dt 



(20) 



Brunone et al.fSa] identify a with the speed of the pressure wave, while the EIT derivation of eq. (PO]) 
suggests it should be chosen as the mean velocity U itself [s^l • In practice, both ks and a are often regarded 
as a free parameters [H, . 

The signs in eq. pO)) are also a matter of some controversy [s^ HO, \EM ■ Among several proposals we 
mention 12311 



dU 

Tt [z, t] = kBpR { sign [[/] — + a 



dU 



dz 



(21) 



The Brunone constant ks may be derived from eqs. ([T^ and (|16p if we require that CM and lAM agree for 
a 2-independent flow with constant acceleration 



R 



dt W \t\ 



AW -^0.06 (I) 



0.185 



(22) 



2.2. The original momentum transfer model 



Let us consider a stationary flow within a straight pipe of circular section and radius R. Let z be the 
coordinate along the pipe. The flow may be decomposed into mean flow and fluctuations as U'' ~ Uz'p + vlP , 
where z^ is the unit vector in the z direction. U is defined by the condition that A (p) U gives the fiux 
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across the section of the pipe, where A is the area of the cross section and () is the average over the section. 
Henceforth we shall drop the angle brackets on (p) 

In the central region, each scale S is associated to a velocity u[d,U]. The MTM claims that the stress at 
the wall is given by 

T = pUu [S*R, U] (23) 

5*R is a characteristic length which is equal to the Kolniogorov scale for flow in smooth pipes (leading to 
the Blasius Law) or else to the pipe roughness eR for rough pipes at large Reynolds number, leading to the 
Strickler Law [y]. We shall assume the latter in what follows. For the recovery of the Poiseuille regime see 

M- 

The velocity u [6*R, U] which appears in eq. is defined from the mode decomposition of the turbulent 
energy 

/•OO 

u[eR,U]^ = {u^[eR,U]) ^2 dkE[k] (24) 

Of course to apply this formula, we need to know the spectrum of turbulent fluctuations for turbulent pipe 
flow. Gioia and Chakraborty get around this difficulty by assuming that the actual spectrum, at least in the 
relevant short wavelength sector, is the same as for an unbounded homogeneous flow, namely the Kolmogorov 
spectrum E [k] = Cif C^^'^fc"^^^, where Ck ~ 1-5 is the so-called Kolmogorov constant [S^- (Although this 
simple picture must be modified in the dissipative and energy ranges, these corrections are not relevant at 
large Re yno lds numbers and we shall not discuss them explicitly; [6| works with a more realistic spectrum, 
see also [ij]). Even after this drastic simplification, a crucial input is still missing, namely the value of C, 
which is the energy flux feeding the Richardson cascade. ^ is usually written in terms of a turbulent velocity 
scale ut as 

C=f (25) 

For UT, Gioia and Chakraborty assume a linear dependence on the mean velocity 

UT = kU (26) 

It is interesting to contrast eq. (|26p . for example, with the usual expression for the turbulent velocity scale 
in terms of a mixing length L, namely ut ~ LdU/dr [l|, The goal of the mixing length approach is 
to obtain the turbulent velocity scale from local properties of the mean flow. An unbounded homogeneous 
mean flow would excite no turbulence, out of galilean invariance, and so mean flow derivatives must be called 
forth. To the contrary, eq. ([25]) is not intended as a local relationship, but as a global one. In particular, 
the boundary condition at the wall plays an essential, if implicit, role. If we could solve the full problem of 
turbulent pipe flow, we expect to find (if only out of dimensional reasons) something like eq. (j26p at least 
close enough to the pipe axis (Ksl - fssj . Because we cannot truly derive eq. (|26|) . we cannot estimate the value 
of the parameter k. However, if we leave it as a free parameter we may compute 

u [eR, Uf 3CKK^e^^^U^ (27) 
then eq. ([^ reduces to eq. ([1]), provided we identify 



Comparing with eq. ([5]) we obtain k — 
free from ad hoc assumptions. 



/ - 8y3C^«:ei/3 (28) 
0.008. We see that even in its original formulation the MTM is not 
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2.3. Stochastic approach to the MTM 

It is clear that the full Navicr- Stokes equation is too complex for analysis, unless numerically. To make 
progress, we shall substitute the Navier-Stokes equation by a linear stochastic one, devised to give the right 
spectrum of turbulent fluctuations, and therefore enabling us to compute u [ei?, U]. 

Let us begin by Fourier decomposing the fluid velocity 

^e^^-<[t] (29) 

We may think of the Fourier components [t] as the velocity associated to eddies of size at time t. We 
conceive the dynamics of these velocities as the balance of two processes: on one hand they draw energy and 
momentum from larger eddies (and ultimately the mean flow), on the other they are subject to dissipation 
from interaction with smaller eddies [s^, [gO]- We model the first process by a stochastic driving force, the 
second one by an eddy viscosity. In steady homogeneous turbulence both processes balance each other, much 
like in the fluctuation-dissipation theorem of linear nonequilibrium thermodynamics (sij . 
We postulate for the Fourier components a dynamic equation 



d 

81 + ''' 



K[t]=F^[t] (30) 



This means that at any time t the velocity of eddies of size k ^ undergoes a random increment F^dt over 
a lapse dt, which subsequently decays exponentially with a mean lifetime cr^^. The velocity increments are 
Gaussian, and uncorrclated if we look at eddies of different sizes and/or at different times 

iK [t] K' n = (2^)' ^ (k + k') 5 (t - t') A^'iV, (31) 
A^'^ is a projector that enforces incompressibility 

Ar = - ^ (32) 

A representation like this may be derived from the functional approach to turbulence, where the left hand 
side of cq. ([50)) is identified as the inverse retarded propagator, and the self-correlation eq. ([31]) is given 
by a self-energy [52], [g^. We shall be content to propose simple expressions for ak and to reproduce the 
known turbulent spectrum. 

In the inertial range, we expect ak and to depend on the only dimensionful parameter C, from eq. \2b\ . 

Let us assume the linear ansatz cq. (j26p . On dimensional grounds Uk oc {k'^C}^ [63| . An important insight 
from ref. [isj is that to obtain the Virk Asymptote for drag reduction by polymer additives, it is necessary 
to assume that (Tj, is inversely proportional to Reynolds number 

.. = ^(*^<)"' (33) 

where A is a pure number. An important consequence of eq. p3p is that is actually independent of the 
local mean velocity U 

au = \lo {Rkf^ (34) 
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where w is defined in eq. (jl7p . This will entail an important simplification in what follows. 

On the other hand, solving the Langevin eq. ((30)) gives the "fluctuation dissipation theorem" k'^Nk = 
E [k] CTfc and therefore 

Nk = CKOkC"^k-^^l^ = \c^x^^J^ (35) 

Observe that A^fe oc U"^ or else oc U. This is the analog to eq. ([25)1 in the stochastic approach. We shall 
see presently that this relation must be generalized for time-dependent flows. 
This concludes the review of the stationary case. 

3. FROM MTM TO lAM 

We now face the problem of generalizing the above effective linear model to the case where the background 
velocity U is position dependent, though time independent [l^. We cannot reduce this problem to the 



homogeneous one, for example, by considering turbulence in the presence of a mean velocity gradient j64| . 
because, as we already remarked, the relationship of the turbulent fluctuations to mean velocity is not purely 
local. Therefore we shall proceed by developing a natural, but to some extent arbitrary, generalization of 
the stochastic approach to the MTM. 

The basic idea is that at any time, and within the volume of fluid between two cross sections at positions 
z and z + eR, the flow remains traslation invariant enough that we may still introduce a Fourier trans- 
form, though the Fourier amplitudes now carry a "slow" variable indicating where the flow is being Fourier 
analyzed. In formulae 

n'M= e^-ul[t;z] (36) 

jRk>e-i (2tt) 

We could make this definition precise by introducing centroid and relative variables (65| . but an intuitive 
idea is sufficient for present needs. See also [38, 62]. Please observe that we shall never need to compute 
other than instantaneous velocity correlations between two points at a distance less than eR apart. 

To obtain a dynamics for the Fourier amplitudes in eq. p6p we assume that, besides the dissipation and 
noise already accounted for in eq. ([5(1)) . turbulence is advected by the mean flow [4l| - |43| . Therefore we write 



- + t/(z,i)-+a. 



K[t;z]^F^[t;z] (37) 



When the Langevin equation is systematically derived from the Schwinger-Dyson equations, it is seen that 
the left hand side of eq. ([37]) is the inverse of the causal propagator for linearized fluctuations in the turbulent 
flow. The evolution of these fluctuations is determined mostly by their interaction with the turbulent eddies 
and it is robust with respect to changes in the external conditions. Therefore it makes sense to assume that 
the mean lifetime of a fluctuation is the same as in a fiducial steady flow. Then ak is still given by eq. p4p , 
and it is position and time- independent. The solution to eq. ([37]) is 



ui[t;z]= f dt'e-'^''i'-'')FP[t';^[z,t;t']] (38) 

J — oo 

where ^ obeys 



^[z,t;t']=0 (39) 
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with boundary condition 



(40) 



To see the meaning of ^, observe that we could define a function zq (t) by holding t' and ^ [zq (t) , t; t'] = 
constant, and then we get dzo/dt = U [zq (t) ,t\ and zq [t') ~ so ^ is a Lagrangian coordinate for the 
particles, with respect to the mean flow. 

Wc now generalize the MTM ansatz eq. ((23)) to inhomogeneous situations as 



T (z, t) = p (z, t) U (z, t) u (z, t) 



(41) 



where 



u (z, t) 



dk' 



Rk>e-^ {2T:f 7i?fc'>£-i (27r)^ 



(42) 



where x is any point close enough to the cross section through z. To compute the expectation value we need 
to make some assumption regarding the noise self-correlation. We assume the random accelerations [t; z] 
are still Gaussian and uncorrelated at different scales and/or times, and that they scale as the local mean 
flow velocity. This leads to 



(K ^] 2']) - (2^) s (k + k') 5 {t - t') A: 



U it,z)U {t,z')k~ 



(43) 



where we have build in the requirement that it reduces to eq. (|43p if the mean velocity is homogeneous. 
Using eqs. (j43l), dM]), ([Ml) and ([Ml) we get 



u (z, t) 



2 ff 



dx 

X>1 ^ 



— f ndt'e-'<'-''yU^[^[z,t;t'],t'] 



(44) 



where 



= eRk 



(45) 



and 



(46) 



Observe that 



dx 



X>1 ^ J —oo 



mt' e-f^(*-*')-''' = 1 



(47) 



Therefore, for steady flow we recover eqs. ([T]) and ([5]). In the general case, we have t = Tg + Tt- Therefore, 
ii T = pUu and Tg = (//8) pU'^, then 



1 



(48) 
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Taking the square of this equation and equating to eq. (pi)) we get 



P 2pTs 



1 (I 



f ^ r ndt' e-^'i'-^'y K [z, t; t'] , t'] - [/2 [z, t]] 

J x>l J -oo 



(49) 



We consider a situation where the mean flow is only weakly inhomogeneous, so the second term in the left 
hand side of eq. (|49l) is negligible, and also time independent. Since z is also ^ we may approximate 



[/^ [C [z, t; t'] , t'] - [z, t] = -{t- t') 2U 



dU d£,[z,t;t'] 



dz dt' 



(50) 



where 



dt' 



[z,t-t'] 



dt 



= U{z,t) 



[z,t;t'] 



dz 



= U{z,t) 



(51) 



P 



16 / ~n~d^ 



(52) 



We recover the original Brunone formula eq. (|20p specialized to time-independent flows, provided we identify 



' 16 / i?0 



(53) 



4. FROM MTM TO CM 



We now consider the complementary problem where the mean velocity is homogeneous in space but time- 
dependent. As in the previous Section, there is no simple way of reducing this problem to the homogeneous 
one. For example, to look at the flow from a comoving, non inertial frame [6^ would be of no avail, because 
it changes the boundary condition at the wall in an essential way. On the other hand, we shall see presently 
that the ansatz eq. (|43p for the noise self-correlation leads to an untenable prediction for transient friction. 
If we believe in the MTM, we must conclude that a time-dependent mean field excites turbulence in an 
essentially different way that a time-independent one. 

Concretely, in the model we have considered so far, the turbulent velocity at a given scale and time is 
the sum of velocity increments occurred at all earlier times and upstream positions. Once produced, each 
elementary increment is exponentially damped while it is advected by the mean flow. 

Suppose wc had a spatially homogeneous accelerating flow. Let us look at the turbulent velocity at some 
time tnow , when the mean flow velocity is Unow ■ If the random velocity increments scale as the instantaneous 
velocity at the time and position of production, as we have assumed so far, then the velocity increments in 
the past have been smaller than for a steady flow with velocity Unow Add the exponential damping to this, 
and the conclusion is that the turbulent velocity, and so the wall friction according to eq. (|4ip . will be less 
for the accelerating flow than for the steady flow with velocity Unow This is contrary to Vardy and Brown 
CM's and most other models in the literature. 

To save the MTM wc must conclude that the random velocity increments in the right hand side of cq. (|37|) 
are sensitive not only to the local mean velocity, but to mean velocity derivatives as well. Of course, lacking 
a full solution for pipe turbulence in time-dependent conditions, our only guidance is phenomenology. On 
the other hand, not to sacrifice too much predictive power, we would like to build a model with no more 
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free parameters than, say, the Vardy-Brown CM. The Vardy-Brown weighting function has two adjustable 
parameters, the A and C constants eqs. (|18p and ([T9|) . Our model already has one parameter over and above 
the steady state model, namely the constant A in cq. (|34p . so we only have room for one new parameter. 

The simplest hypothesis is that the local random velocity increments scale not as the instantaneous 
velocity C/ [z, i] but rather as some effective velocity scale C/ [z, t] which depends also on the local acceleration. 
Repeating the arguments in the previous Section we arrive at 



1 ff 



2/OTs 3 



dx 



x>l 



[t',x\-U' [t] 



(54) 



We impose the requirement that for a weakly time-dependent and space independent flow, eq. (f54| . after 
linearization in Tt and dU/dt, should reduce to eq. ([M]) with a weighting function consistent with Wvb as 
given in eqs. (jlSp and (|16p. When Tt is small, we may write 



where 



Tt 



(55) 



.(1) 



3 VS 

-(^ 

3 \8 



Jx>l J -oo 

J x>l J —oo 



dU 
di/ 



(56) 



.(2) 



P 



1 ff 

3 



dx 



X>1 



ndt' e-"(*-*')-'' 



[t',x]-U^ [t'l 



It is clear that if we choose U = U then Tj = and we get an expression for Tj with the wrong sign, as 
expected. Therefore we must contemplate a more general ansatz. To keep the number of new parameters to 
a minimum, we are satisfied with proposing a simple scale-free form 



(57) 



dU 

U%x]^U[t]+T{Rkr — [t] 

ot 



(58) 



T is a characteristic time to be determined, and it is the single new parameter we shall allow ourselves. After 
linearization we get indeed eq. (|14p with a weighting function 



-X — X 



-2/3 



(59) 



To obtain [t] w t ^/2 when t — > requires a = 1/3. Integrating over time we obtain the Brunone constant 





Re oj 


flT 1" 


a) 


n 


" 4_ 



(60) 



and then, writing x^^"^ = y/y/Tti 



W[t] = { ^]ReRuj sign [U] 



nt Ei [-nt] 



(61) 
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where 



erfc [x] 



dy e~ 



(62) 



Ei [x] 



X y 



(63) 



This weightmg function decays as exp(— f2t). To match the Vardy-Brown weighting function eq. (jl6p we 
ask 

Q^Cu (64) 
where C is defined in eq. ((T^ . This may be used to fix the constant A in eqs. ([55]) and namely 



(65) 



Observe that if we adopt eq. then we move beyond the framework of [T^, because there A is supposed 
not to depend on Reynolds number. It should be observed that ref. [lB| deals with relatively low Reynolds 
numbers, consistent with Blasius' Law in absence of the polymer, while here we are considering very high 
ones to which Strickler's scaling applies. 

Using eqs. ^ and dM]) into §U^i we get 



^"4~ 8^" V2. 



(66) 



For example, ife/2 = 10 ^, then the right hand side of eq. gives 0.15, and becomes smaller by a factor 
of 10 if e/2 = 10~*^. For these values of flT/e^/^ the MTM weighting function eq. (|6T|) changes sign for 
relatively low values of fit. Nevertheless, the scale-free ansatz eq. (|58p captures the short time behavior and 
the overall integral of the weighting function, so it may be trusted to give right results both for strongly 
accelera ting flows and slowly varying flows. In any case, all CM models have only a finite time duration of 
validity |67H69j , so we cannot truly say that such a sign reversal is ruled out by experiment [t^ . 

For larger values of QT/e^^^ the sign reversal occurs too late to be of any relevance. This follows from the 
asymptotic expansions [71[ 



erfc I 



(67) 



Ei [-x] w - — 

X 

After the identifications eqs. ([501) and the Vardy and Brown weighting function eq. ((T5)) reads 



WvB [t] = ( J ) Re Roj 



7^ 



(69) 



In fig. 1 we show a plot of both weighting functions as a function of ilt, divided by (//4)Re Roj, for 
ilT/e^/^ = .4. Fig. 2 is a close up of the short t behavior, for which we have chosen a log-log scale. 
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FIG. 1; [Color online] Plot of the Weighting function W (eq. (|6ip Ufull line) and the Vardy-Brown weighting function 
WvB (eq. Uni) (dashed), both divided by (//4) Re Riu, as functions of ilt 




FIG. 2: [Color online] Log-Log plot of the Weighting function W (eq. (|61(l )(full line) and the Vardy-Brown weighting 
function Wvb (eq. ([16])) (dashed), both divided by (//4) Re Rlo, as functions of Qt 



5. FINAL REMARKS 



In this paper we have shown a simple generalization of the momentum transfer model to time-independent 
inhomogeneous flows, which yields a dynamic friction similar to the Brunone model [36| . 

Matching the Vardy and Brown model [s^ [sjj proves to be a harder task. In this case we must appeal to 
a phenomenological approach to obtain a minimally acceptable MTM model 

Although we have made several ad hoc choices in order to achieve this match, we believe there are at least 
two senses in which these results are meaningful. 

First, if it is conceded that a weakly time dependent mean flow produces a spectrum of turbulent fluctu- 
ations as depicted by the stochastic equation eq. (|37)) with the noise self correlation given by eqs. (|43|) and 
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(|58p . then this may be used as a benchmark for more fundamental approaches such as those in [65f and [6^ . 

Second, while the model is built to match Vardy and Brown's weighting function for a weakly time- 
dependent flow, it is clearly superior to it when space dependence becomes an issue, since it has the Brunone 
model built in. This synthesis of the lAM and CM expressions for unsteady friction is a legitimate prediction 
of the momentum transfer model, which it ought to be possible to contrast against experimental data. 

We continue our research in these two promising directions. 
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